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1 Introduction 

Recently, the high-temperature electroweak phase transition has been the subject of 
active research, due to its possible effect on the baryon number of the Universe. The 
standard tool for perturbative investigations of the equilibrium properties of this phase 
transition is the effective potential for the Higgs field, V{!f). From the effective poten- 
tial, one can calculate quantities like the order, the critical temperature, and the latent 
heat of the phase transition. In particular, the pressure is given by minus the value 
of the effective potential at the minimum. Presently, the most complete expressions 
for V{ip) are given by the 2-loop calculations in refs. [0, H, 0- 

When calculating the effective potential, one of the problems one is faced with is 
gauge dependence. Indeed, the source term appearing in the generating functional 

exp{-W[J]) = fV(f)VAexp{-S[(j),A] - (J^$ + $V)/2) (1) 

is not gauge- invariant, and hence W[J] may be gauge-dependent. The effective po- 
tential V{ip), which is obtained from the Legendre transformation of H^[J], then also 
depends on the gauge condition [^]. In particular, the location of the minimum of the 
effective potential is gauge- dependent, which indicates that ip is not a physical observ- 
able. However, the value of the effective potential at a stationary point is obtained by 
putting the source J to zero, and consequently one should obtain a gauge-independent 
result for the pressure @, . It is the main purpose of this paper to check explicitly to 
2-loop order (^^) that a gauge- independent result is, indeed, obtained. 

There are methods of extending the definition of the effective potential so that it is 
gauge-independent even away from the minima (see, e.g., ref. [0). A concrete example, 
relevant for the EW phase transition, was given in ref. p]. There the external source 
was coupled to the composite operator $^$: 

exp(-i^[J]) = fV(j)VAexp{-S[(f),A] - 2J$^$) . (2) 

The source term 2J$^$, and hence both W[J] and the corresponding effective po- 
tential V{a), are manifestly gauge-independent. The pressure is again obtained from 
the value of V{a) at its minimum, or, equivalently, from iy[0]. It is checked below to 
order h'^ that V{a) and the conventional gauge-dependent effective potential V{(p) give 
the same result for the pressure. 

The paper is organized as follows. In Sec. ^Method of calculationsection.2| , some 
technical details of the calculation of the 2-loop effective potential in a general co- 
variant gauge are discussed. In Sec. |3Equation of statesection"3|, it is shown that 



the value of the effective potential at the minimum is gauge-independent to order h . 
In Sec. |4Gauge-independent ehective potentialsection.4| , it is shown that the effective 



potentials resulting from eqs. |T] and ^ yield the same values for the pressure, when cal- 
culated consistently in powers of h. The calculation of some other physical quantities 



than pressure is discussed in Sec. |50ther physical quantities in powers of ^section.5| . In 



Sec. |6Convergence in the broken phasesection.q , the convergence of perturbation theory 
is studied in the broken phase, utihzing the gauge- dependent effective potential. The 
conclusions are in Sec. |7Conclusionssection. 7| . The explicit form of the 2-loop effective 
potential is represented in the Appendix. Sections ^Method of calculationsection.2 
and ^Convergence in the broken phasesection."^ are an extension of ref . ^ . 



2 Method of calculation 



In the high-temperature limit, the essential long- wavelength properties of the Stan- 
dard Model can be described by an effective super-renormalizable 3-dimensional field 
theory IJIO, 11|. This dimensional reduction is accomplished by integrating out all the 
degrees of freedom corresponding to the momentum scale p ~ T in the original 4d the- 
ory. At the 2-loop level, it has been explicitly verified that the 3d theory produces 
the high-temperature limit of the 4d effective potential for the Higgs field. Further- 
more, the parameters of the effective 3d theory should be independent of the 4d gauge 
fixing condition M, |T2l. The relevant 3d action is 



S 



(TX 



F^^F^^ + (A$)^(A$) + m^$t$ + A3($t$)' + h^^^A^Al 



4r "■! 
1 



+ -{D,AIY + -m'j,AlAl + -M^IK) 



(3) 



where F^, 
''J 



{di - 2^3rM^/2)<|), and A^S = diA^o + 



diA] - djA^ + gae^'^'A'lA'j, Di<l> 
g^e^'^^A'^AQ. The r":s are the Pauli matrices. All the fields have the dimension GeV^/^, 
and A3, g^, /13 and A^ have the dimension GeV. 

The theory in eq. ^ is a simplified theory in the sense that the U(l)-subgroup 
has been neglected. All the fermions can be included, but they appear only in the 
parameters of the effective theory. The relations of the parameters of the 3d theory to 
those of the Standard Model are given in refs. [§, |13|, but we do not presently need 
these relations. In some cases, the 3d theory could be further simplified by integrating 
out also the Ag-field P, Q, resulting in a 3d SU(2)-Higgs model. For generality, we 
keep the Ag-field in the action. In the following, we calculate the effective potential in 
the theory defined by eq. ^. Due to the change in the dimension of the space in the 
reduction step, the value of this effective potential at the minimum is then minus the 
pressure divided by the temperature. 

The Lagrangian in eq. ^ is gauge-invariant in r-independent gauge transforma- 
tions, and gauge fixing and compensating terms have to be added for perturbative 
calculations. We choose these as 



S, 



d^x 



1 
2^ 



\f> „c 



(diA^y + dic'^diC^ + gse'^'^dic'^A^d 



(4) 



The gauge parameter ^ is a renormalized version of the corresponding 4d gauge fixing 
parameter. The covariant gauge condition of eq. ^ has both advantages and disadvan- 
tages compared with other usual gauges, hke the i?^-gauge. On one hand, it is clearly 
an allowed choice of gauge, whereas with the i?^-gauge, one is choosing a different gauge 
for each different value of (p in the effective potential, which might not be allowed. On 
the other hand, the covariant gauge generates extra IR-divergences, as is seen below. 

To calculate the effective potential V{lp), one writes $ = [03+Z04, y9+0i+z02]^/v^ 
in the action S + S^ and neglects terms linear in quantum fields ^. This defines a new 



theory with the masses m\ = ml + SAsy^^ 



mr, 



m\ + Aav^^, m\ = mjj + h^if'^, and 



rrin 



g^ifi /A. The propagators are 
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k"^ + m\ 
k"^ + m?T 
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-k)A]{k) 
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5, 



kikj / k ^kikj 



k"^ + ml 



M-k)Mk) 
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Mk)AK-k) 



k"^ + iTij, 
k"^ + S,m\ 



k'^ik'^ + m\) + ^m^m^ ' 
-03(A:)A'(-^) = -Uk)A]{-k) 



k"^ k"^ {k'^ + ml) + ^rrij^ml 
G = 2,3,4 

i^TTlTk, 



(5) 



A;2(A;2 _|_ ^2^ _|_ ^m'^nri'l 



If ml is negative, a small imaginary part has to be added to it to define the loop 
integrals. One then calculates all the one-particle-irreducible vacuum diagrams of the 
theory to a desired order in the loop expansion. Non-vanishing 2-loop contributions 
arise from the diagrams of Fig. 1. The method of calculation is to write 



1 



1 



k^ {k^ + ml) + ^m^ml mliR^ — R 



1 



1 



A;2 + mlRl k'^ + mlRl 



(6) 



where R\. = 1/2 ± \/l/4 — ^(771^/7712)^, and to use standard Landau-gauge values of 
integrals j^, with dimensional regularization. The 2-loop effective potential is presented 
in the Appendix. 

Next, the effective potential has to be renormalized. Many of the individual 2-loop 
diagrams contributing to V{ip) include gauge- dependent divergent pieces, but all these 
cancel, leaving the gauge-independent divergence 



h'^ ji 



-4e 
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^3 + 9A3^3-12A2 + 12/13^2 



Qhl 



?>glml + Qglm]^ \ . (7) 



The piece multiplying Lp^ is removed by mass renormalization, and the coupling con- 
stants are RG-invariant. Note that if the vacuum terms are renormalized by just 



removing the 1/e-pieces, the value of the effective potential at the minimum be- 
comes /i-dependent. The yj-dependent part of V{ip) is /^-independent, since the /i- 
dependence of the renormalized mass squared rnKfi) cancels the log(/x)-terms of the 
2-loop graphs. The renormalized effective potential is obtained from the expressions 
in the Appendix by replacing m'^ with 7n|(/x), and by ignoring the 1/e-piece in the 
function H{ma,rnb,mc). 

3 Equation of state 



To calculate the value of V{(f) at the minimum in powers of h, we write V{(f) as 
V = Vq + hVi + ^^¥2, with the classical part Vq = m|(^^/2 + X^ip^/A. The location of 
the minimum is determined from 



d 

dip 



Vo{^)\ 



(p=ipo+hipi+h ip2 



+ ^^l(^)Uo+n^,+^' ^2(^)1 



ip=ipO 



The ^°-term of this equation reads Vo(v5o) = fo^ = 0, yielding the two solutions 
(po = and ip^ = —m\/\-i. Let us assume that m\ < 0, and inspect first the broken 
minimum. From eq. ^, the corrections to ipo are 

1 {v;r 



Vl 






^2 = - 



1 d 
V^'Jp 



Vo- 



2 K," 



(9) 



where each expression is evaluated at ipQ. The value of the 3d effective potential at the 
minimum is then 



P{T) 
T 



K)(v5o) + ^ 



Vi + <^iK,' 



+ K' 



IP = IPQ 



2 K" °. 



(10) 



¥'=¥'0 



This should be gauge- independent. 

At 990 ) K) ~ t)y definition, so that the terms multiplying (pi and (p2 in eq. |T^ 
should be put to zero. However, these terms are kept for the moment, for the following 
reason. In the limit m^ — ^ the Goldstone mode propagator is 



1 im% 



A;2 + 



k^ 



:iii 



Because this is IR-divergent inside loop integrals, the 2-loop effective potential diverges 
in the limit m\ —>■ 0. The divergences from the individual diagrams which are of the 



form ^5/2^^\ ^^/^ 



771. 



-1/2 



or ^2log(m2), cancel, but the divergent piece 
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remains. These problems show up even at the 1-loop level: the ^-dependent part of Vi 
is finite but non-analytic, 

V;«) oc -e^/^m^/' . (13) 



Due to the equation dm^jd^ = rm2~'^\^^, it is seen from eqs. ^, ^ and |T3| that the 
correction y^i diverges as —^^/^m2 , and Lp2 might diverge as ^^/"^m'^^ . However, it 
can easily be seen from eq. |2^ that the term 

-i^ (14) 

exactly cancels the divergent piece of V2, shown in eq. [l^. Consequently, it turns 
out that Lp2 only diverges as — ^m^\ so that being multiplied by Vq oc m\ in eq. |iy, 
it does not contribute to the value of the effective potential at the minimum, where 
m\ = 0. In other words, the terms f/^iV^' and ip2Vf^ can safely be put to zero, when 
the IR-divergences are "regularized" by handling the classical minimum as the limit 
1712 -^ 0^- Let us mention that in the Landau-gauge (^ = 0), (fi and (f2 are finite. 

In addition to cancelling all the divergences, the term in eq. |l^ also cancels the 



finite gauge-dependent piece of V2 in eq. |T0[ As mentioned, the gauge dependence of 
Vi is of the form ^I'^m2 near Lp^^ so that this vanishes, too. Thus, we have verified 
explicitly that the pressure at the broken minimum is gauge-independent to order T? . 
At the symmetric minimum </? = 0, to 2-loop order there are no IR-divergences 
in V2, and the effective potential is gauge- independent in itself. The "correction term" 



in eq. |T4| vanishes. If m\ < 0, the effective potential becomes complex, in accordance 
with the fact that the symmetric minimum is, at the tree level, unstable for m"^ < 0. 

In addition to C,, all the physical quantities should be independent of the renor- 
malization scale /x. As mentioned after eq. |^, in dimensional regularization the di- 
mensionally reduced effective theory produces a /^-dependent unphysical vacuum term 
to the 3d effective potential, corresponding to a zero-temperature entropy. This term 
is actually removed by a counterterm produced by the dimensional reduction step, 
but on the other hand there is an undetermined vacuum energy density in the full 4d 
theory. Apart from the vacuum term, eq. ^ is /x-independent to order h , since the 
yj-dependent part of V{ip) is so. In particular, the difference between the pressures of 
the symmetric and broken phases, determining the critical temperature and the latent 
heat of the phase transition, is independent of both /i and ^ to order h'^, as it should. 

4 Gauge-independent effective potential 

Let us calculate the gauge-independent generating functional defined in eq. ||, and the 
value of the corresponding effective potential at the minimum, in powers of h inside 
the 3d theory. Since we are only interested in the effective potential, J may be chosen 
as a constant. Defining w{J) = W[J]/V and dw/dJ = a, the value of the effective 
potential V{cr) = w — a J at the minimum is formally 

'^WLvV^.o = -(»)■ (15) 



To see how this equation arises in powers of h, write w{J) = wq + hwi + h'^W2. For 
a Legendre transformation, the source J = Jq + hJi + h^J2 has to be solved from the 
equation dw/dJ = a. With Jq determined from Wo(Jo) = <7) the effective potential is 
then 

V{(y) = Wo{Jo) - aJo + hwi{Jo) + ^^ 



W2 



2 w'^ 



(16) 

J=Jo 



Calculating the value of V{cr) at the minimum as in eqs. ^^, finally gives the result 
l/(min) = wo{0) + hwi{0) + h'^W2{0) as in eq. |^. Below, we see exphcitly that wo{0) + 
hwi{0) +h'^W2{0) is given by eq. [1^, so that as far as the equation of state is concerned, 
the conventional effective potential V{ip) and the gauge-independent effective potential 
V{a) should lead to the same physics. 

To calculate wlJ) [and, in particular, ti'(O)], we need to perform the path integral 

fV(f)VAexp{-S[(t),A]) (17) 

using the loop expansion, and then to substitute ml —* ml + 2 J. We first discuss the 
broken minimum. Let us make the change of variables $ = (0, ipo)'^ / \^ + $' in eq. |T^, 
where ipo = J—m^l/X-^ is the location of the classical broken minimum. Since this 
is a stationary point of the Lagrangian, all terms linear in quantum fields disappear 
from the action (apart from counterterms contributing to higher order than h ). This 
means that we have exactly the same Lagrangian as in the case of calculating the 
effective potential V{ip) at ip = ipo. However, the set of diagrams is different, since 
in the calculation of the effective potential, only one-particle-irreducible graphs are 
included ^ . Now we have to include all the connected graphs, which means that one- 
particle-irreducible graphs have to be supplemented by graphs of the type shown in 
Fig. 2. These graphs can be easily calculated in the same gauge which was used before. 
The result is simple: these graphs contribute the amount — (V"/)^/(2Vq'). Hence, w{0) 



equals the gauge-independent value of eq. |I0|, as promised. 



When constructing the gauge-independent effective potential in the symmetric 
phase, there is the complication that the tree-level part of the generating function 
w{J) vanishes, Wo{J) = 0. Therefore, the calculations after eq. [l^ have to be modified. 
Nevertheless, it is still true that w{0) as calculated from the loop expansion equals the 
value of the conventional effective potential at ip = 0. 

5 Other physical quantities in powers of h 

Having constructed an expression for the value of the 3d effective potential at its 
minima, which is independent of ^ and independent of /i to order h , we could express 
the 3d parameters in terms of the 4d parameters and temperature, and study the 
properties of the EW phase transition. Of course, it is expected that this investigation 
breaks down at some order, due to the IR-divergences of finite-temperature field theory. 

6 



Let us see, how far we can go without problems. From the difference Ap{T) of the 
pressures of the symmetric and broken phases, we can in principle calculate the critical 
temperature T^ and the latent heat L to order fi^. The former is the solution of the 
equation Ap{Tc) = 0, and the latter is given by L = TcAp'{Tc). Unfortunately, the 
zeroth order term Tq of the critical temperature Tc = Tq + hTi + h?T2 is given as the 
solution of the equation 

m^(To,/i) = 0, (18) 

and the vanishing of m\ leads to singularities at high enough orders. The first order 
term Ti is still finite, and given by 



Ti 



Api(To 
Ap[,(To 



A-K 



h^niD 



It 



(19) 



T=To 



To leading order in the coupling constants, this agrees with the result in ref. |T5 



However, the second order term T2, which includes a finite piece cancelling the fi- 



dependence of Tq arising from eq. 0, also includes a logarithmic divergence T2 oc 
log(yu/r?7,3). Hence, Tc is not calculable to order h . As to the latent heat, both the 
zeroth order term Lq = TqAp'q[To) and the first order term 



Li = -Api(To)+To Ap'i- 






(20) 



T=To 



vanish. For Lq this is natural, but for Li it is not, since the 1-loop effective poten- 
tial V{ip) in the Landau gauge has a barrier between the symmetric and broken minima, 
suggesting a first-order phase transition. The second order term L2 is finite (after can- 
cellations of various divergences, proportional to m^'^, m^^ and logms), but it includes 
an imaginary part: 



rp2 

167r2 




A3 
2 



4A3 V16^^ ^' 




(21) 



T=To 



The real part of L2 behaves qualitatively like the (7^, A^ -curve in Fig. 6 of ref. [0, but 
the numerical value is much smaller. We conclude that little constructive information 
concerning the actual EW phase transition can be obtained from the present approach. 
Even if the physical properties of the EW phase transition cannot be reliably 
calculated in powers of h, it might be hoped that quantities related solely to the broken 
phase, away from the critical temperature, could be calculated. The problem is that 
there may be //-dependent vacuum parts in these quantities as in the value of pressure. 
However, when one is comparing perturbation theory with lattice calculations, fi can 
be fixed |l^, and the vacuum parts are not a problem. To illustrate the convergence of 



these calculations, we note that for mn = 80 GeV, /i = 100 GeV and T = 60 GeV, the 
^°, n^ and h'^ contributions to ($"f$) are [25.47, 3.76 and 0.39] x 10^ GeV^ (for these 
parameters, T^ ~ 170 GeV). These values seem to indicate reasonable convergence. In 
the next Section, the convergence of perturbation theory in the broken phase is studied 
from a different point of view. 

6 Convergence in the broken phase 

Above, we were careful to calculate quantities which are gauge-independent to each 
order in h, and we could just hope that this expansion converges. One could also take a 
different approach, trying to extremize convergence but not being too careful about the 
gauge dependence. The possible gauge dependence of the physical observables obtained 
in this way could then be used as a measure of the convergence of the expansion. 

Good convergence of perturbation theory requires that the coupling constants, and 
the higher order contributions including logarithms of the renormalization scale n, are 
small. The latter requirement can be satisfied by making a renormalization group (RG) 
improvement to the naive calculation. In practice, this can be implemented by choosing 
II to correspond to a typical mass scale appearing in the propagators. When calculating 
the effective potential in a large range of (p, the masses in the propagators depend on 
the shifted field ip, and hence fi should also be chosen to depend on ip. When this is 
consistently implemented, one arrives at the RG- improved effective potential |^, |T^. 

In ref. 0, the physical observable dp /dm? = — ($^$) was chosen as the indicator of 
the convergence of the RG-improved loop expansion. The value of ($'''$) was calculated 
numerically from the RG-improved 2-loop effective potential. Due to RG- improvement, 
and the fact that one is calculating the location of the broken minimum exactly instead 
of using eq. |^, the result includes a certain subset of contributions from higher powers 
of h. Since {^^^) is gauge-independent in the full theory, the gauge dependence of 
the mixed-order result tells something about the order of magnitude of the missing 



terms. It is seen in Fig. 3 of ref. that w2($t$) depends much less on the gauge 
fixing parameter than, for instance, the location of the broken minimum, which is not 
gauge-independent. To give a numerical illustration, for ^ = and the parameters 
cited at the end of Sec. |50ther physical quantities in powers of /isection.5| , the value of 



($'''$) is 29.64 X 10^ GeV^, and the uncertainty due to gauge dependence is of order 
1 percent. This gives support to the argument that the RG-improved loop expansion 
converges well deep in the broken phase. 

7 Conclusions 

By calculating consistently in powers of h, we have been able to derive a gauge- 
independent expression for the pressure of EW matter at high temperature. Unfortu- 



nately, the physical properties of the EW phase transition cannot be rehably calculated. 
However, quantities related solely to the broken phase are calculable in powers of h, 
and the convergence of perturbation theory should be reasonably good. 
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Appendix 

We present here the 2-loop effective potential of the theory defined by eqs. ^ and ^. 
The tree-level part is 

V,{v) = \mlv^ + \W ■ (22) 

With the functions R± defined in eq. ^, the 1-loop contribution to the effective potential 

'' 1 

Vii^) = -T7^[6^T + 3mi + m\ + ?,ml{Rl + RI)] . (23) 

In dimensional regularization, there are no divergences in Vi{ip). To present the 2-loop 
contribution, we use the function 



H{ma,mb,mc 



167r2 



+ log' '^ ^ ■ 



4e \ma + nib + mcJ 2 



(24) 



arising from the sunset diagrams in Fig. 1. For brevity, we also use the functions 
I{ma,mf,,mc), Dvvs{M,M,m) and Dssv{'>TT'a,^TT'b,M) defined in ref. [Q. These func- 
tions are shorthands for certain combinations of masses and logarithms, and the func- 
tion H{ma,mb,mc) appears linearly in them. Four more functions, denoted by A, B, 
C and V, are defined at the end of this Appendix. A common factor l/lGvr^, and the 
factor /i"^"^ multiplying H{ma, mi,, rric), have been omitted from all the formulas below; 
accordingly, the factor l/lGvr^ in eq. ^ should also be left out. The divergent part of 
V2{(p) is obtained from the coefficient of the function H{ma,mb,mc), and is displayed 
in eq. |^. To get the renormalized expression for V2{(p), the term 1/e is to be omitted 
from eq. ^ and m| is to be replaced by mKfi). Let us note that due to the identities 
R"^ + R^ = 1 and R^R^ = ^{mT/rn2)'^, it might be possible to simplify some of the 
formulas below. We denote m2-R± by m±. Finally, for pure SU(2)-Higgs theory without 
the ^Q-field, the graphs (g3), (e2), (g4), (b) and (el) are left out. 
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mi 2 + 
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— m^+ 217111712^ 52" + ^"^2 
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^■' ' -- mr mi + 3m2-r± ^2 + 
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/(o. 0. ,„o - /(o. ,„,. ,„,) - /(o. ™_. ™o + /(,„_. ™,. ™,)1 + (+ " -)■ 

i/(mi,mL,mL) 

i/(mi,mi,mi) 
1 



(i?i - i?^)2 



R^H{m^, m_|_, mi) — R_R_^H {171^,711^, mi) + (+ ^^ — ) 



(el) 
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Fig. 1: The non-vanishing 2-loop graphs. Dashed line is the scalar propagator, wiggly 
line is the yl"-propagator, solid line is the ylo'Propagator, and double line is the ghost 
propagator. The notations are the same as in ref. [3]. 




Fig. 2: The general structure of the connected reducible diagrams that are not included 
in the calculation of the effective potential, but must be included, if one is calculating 
the partition function in the loop expansion without using the effective potential. Only 
Higgs particles can propagate on the straight line. 



